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Abstract — In this paper, we provide a throughput analysis 
of the IEEE 802.11 protocol at the data link layer in non- 
saturated traffic conditions taking into account the impact of 
both transmission channel and capture effects in Rayleigh fading 
environment. The impact of both non-ideal channel and capture 
become important in terms of the actual observed throughput in 
typical network conditions whereby traffic is mainly unsaturated, 
especially in an environment of high interference. 

We extend the multi-dimensional Markovian state transition 
model characterizing the behavior at the MAC layer by including 
transmission states that account for packet transmission failures 
due to errors caused by propagation through the channel, along 
with a state characterizing the system when there are no packets 
to be transmitted in the buffer of a station. Finally, we derive a 
linear model of the throughput along with its interval of validity. 

Simulation results closely match the theoretical derivations 
confirming the effectiveness of the proposed model. 

Index Terms — Capture, DCF, Distributed Coordination Func- 
tion, fading, IEEE 802.11, MAC, Rayleigh, rate adaptation, 
saturation, throughput, unsaturated, non-saturated. 



I. Introduction 

WIRELESS Local Area Networks (WLANs) using the 
IEEE802.il series of standards have experienced an 
exponential growth in the recent past [1-23]. The Medium 
Access Control (MAC) layer of many wireless protocols 
resemble that of IEEE802.il. Hence, while we focus on this 
protocol, it is evident that the results easily extend to other 
■ protocols with similar MAC layer operation. 

The IEEE802.il MAC presents two options (T], namely 
'the Distributed Coordination Function (DCF) and the Point 
Coordination Function (PCF). PCF is generally a complex 
access method that can be implemented in an infrastructure 
.network. DCF is similar to Carrier Sense Multiple Access 
with Collision Avoidance (CSMA/CA) and is the focus of this 
paper. 

With this background, let us provide a quick survey of the 
recent literature related to the problem addressed here. This 
survey is by no means exhaustive and is meant to simply 
provide a sampling of the literature in this fertile area. We 
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invite the interested readers to refer to the references we 
provide in the following and references therein. 

The most relevant works to what is presented here are lfT3l . 
llT4l . ifTSl . In IIT3I the author provided an analysis of the 
saturation throughput of the basic 802.11 protocol assuming 
a two dimensional Markov model at the MAC layer, while 
in [14 1 the authors extended the underlying model in order 
to consider unsaturated traffic conditions by introducing a 
new idle state, not present in the original Bianchi's model, 
accounting for the case in which the station buffer is empty, 
after a successful completion of a packet transmission. In 
the modified model, however, a packet is discarded after m 
backoff stages, while in Bianchi's model, the station keeps 
iterating in the m-th backoff stage until the packet gets 
successfully transmitted. In [15], authors propose a novel 
Markov model for the DCF of IEEE 802.11, using a IEEE 
802.11a PHY, in a scenario with various stations contending 
for the channel and transmitting with different transmission 
rates. An admission control mechanism is also proposed for 
maximizing the throughput while guaranteeing fairness to the 
involved transmitting stations. 

In [16 1, the authors extend the work of Bianchi to multiple 
queues with different contention characteristics in the 802.1 le 
variant of the standard with provisions of QoS. In ifTTl . 
the authors present an analytical model, in which most new 
features of the Enhanced Distributed Channel Access (EDCA) 
in 802. lie such as virtual collision, different arbitration in- 
terframe spaces (AIFS), and different contention windows 
are taken into account. Based on the model, the throughput 
performance of differentiated service traffic is analyzed and a 
recursive method capable of calculating the mean access delay 
is presented. Both articles referenced assume the transmission 
channel to be ideal. 

In [18 1, the authors look at the impact of channel induced 
errors and the received SNR on the achievable throughput in 
a system with rate adaptation whereby the transmission rate 
of the terminal is adapted based on either direct or indirect 
measurements of the link quality. In [T9l, the authors deal with 
the extension of Bianchi's Markov model in order to account 
for channel errors. In ll20l . the authors investigate the satura- 
tion throughput in both congested and error-prone Gaussian 
channel, by proposing a simple and accurate analytical model 
of the behaviour of the DCF at the MAC layer. The PHY-layer 
is based on the parameters of the IEEE 802.11a protocol. 

Paper [21] proposes an extension of the Bianchi's model 
considering a new state for each backoff stage accounting for 



the absence of new packets to be transmitted, i.e., in unloaded 
traffic conditions. 

In real networks, traffic is mostly unsaturated, so it is important 
to derive a model accounting for practical network operations. 
In this paper, we extend the previous works on the subject 
by looking at all the three issues outlined before together, 
namely real channel conditions, unsaturated traffic, and capture 
effects. Our assumptions are essentially similar to those of 
Bianchi fT3l with the exception that we do assume the 
presence of both channel induced errors and capture effects 
due to the transmission over a Rayleigh fading channel. As 
a reference standard, we use network parameters belong- 
ing to the IEEE802.11b protocol, even though the proposed 
mathematical models hold for any flavor of the IEEE802.il 
family or other wireless protocols with similar MAC layer 
functionality. We also demonstrate that the behavior of the 
throughput versus the packet rate. A, is a linear function of 
A up to a critical packet rate, Ac, above which throughput 
enters saturated conditions. Furthermore, in the linear region, 
the slope of the throughput depends only on the average packet 
length and on the number of contending stations. 

Paper outline is as follows. After a brief review of the 
functionalities of the contention window procedure at MAC 
layer, section II extends the Markov model initially proposed 
by Bianchi, presenting modifications that account for transmis- 
sion errors and capture effects over Rayleigh fading channels 
employing the 2-way handshaking technique in unsaturated 
traffic conditions. Section III provides an expression for the 
unsaturated throughput of the link. In section IV we present 
simulation results where typical MAC layer parameters for 
IEEE802.11b are used to obtain throughput values as a func- 
tion of various system level parameters, capture probability, 
and SNR under typical traffic conditions. Section V derives 
a linear model of the throughput along with its interval of 
validity. Finally, Section VI is devoted to conclusions. 

II. Development of the Markov Model 

In this section, we present the basic rationales of the 
proposed bi-dimensional Markov model useful for evaluating 
the throughput of the DCF under unsaturated traffic conditions. 
We make the following assumptions: the number of terminals 
is finite, channel is prone to errors, capture in Rayleigh fading 
transmission scenario can occur, and packet transmission is 
based on the 2-way handshaking access mechanism. For 
conciseness, we will limit our presentation to the ideas needed 
for developing the proposed model. The interested readers 
can refer to lUll for many details on the operating 
functionalities of the DCF. 

A. Markovian Model Characterizing the MAC Layer under 
unsaturated traffic conditions, Real Transmission Channel and 
Capture Effects 

In lfT3l . an analytical model is presented for the computation 
of flie throughput of a WLAN using the IEEE 802.11 DCF 
under ideal channel conditions. By virtue of the strategy 
employed for reducing the collision probability of the packets 
transmitted from the stations attempting to access the channel 



simultaneously, a random process h{t) is used to represent 
the backoff counter of a given station. Backoff counter is 
decremented at the start of every idle backoff slot and when 
it reaches zero, the station transmits and a new value for h{t) 
is set. 

The value of b{t) after each transmission depends on the size 
of the contention window from which it is drawn. Therefore 
it depends on the stations' transmission history, rendering it a 
non-Markovian process. To overcome this problem and get to 
the definition of a Markovian process, a second process s{t) is 
defined representing the size of the contention window from 
which h{t) is drawn, [W^ = i = s{t)). 

We recall that a backoff time counter is initialized depending 
on the number of failed transmissions for the transmitted 
packet. It is chosen in the range [0, W^ — 1] following a uniform 
distribution, where Wi is the contention window at the backoff 
stage i. At the first transmission attempt (i.e., for i = 0), 
the contention window size is set equal to a minimum value 
Wq — W, and the process s{t) takes on the value s{t) = i = 0. 

The backoff stage i is incremented in unitary steps after 
each unsuccessful transmission up to the maximum value m, 
while the contention window is doubled at each stage up to 
the maximum value CWmax — 2™W^. 

The backoff time counter is decremented as long as the 
channel is sensed idle and stopped when a transmission is 
detected. The station transmits when the backoff time counter 
reaches zero. 

A two-dimensional Markov process {s{t), b{t)) can now be 
defined, based on two assertions: 

1) the probability r that a station will attempt transmission 
in a generic time slot is constant across all time slots; 

2) the probability Pcoi that any transmission experiences a 
collision is constant and independent of the number of 
collisions already suffered. 

Bianchi's model relies on the following fundamental assump- 
tions: 1) the mobile stations always have something to transmit 
(i.e., the saturation condition), 2) there are no hidden terminals 
and there is no capture effect (i.e., a terminal which perceives 
a higher signal-to-noise ratio (SNR) relative to other terminals 
capturing the channel ll22l . Il23l . Il24l and limiting access to 
other terminals, similar to the near-far problem in cellular net- 
works), 3) at each transmission attempt and regardless of the 
number of retransmissions suffered, each packet collides with 
constant and independent probability, and 4) the transmission 
channel is ideal and packet errors are only due to collisions. 
Clearly, the first, the second and fourth assumptions are not 
valid in any real setting, specially when there is mobility and 
when the transmission channel suffers from fading effects. 

The main aim of this section is to propose an effective 
modification of the bi-dimensional Markov process proposed 
in 1,13] in order to account for unsaturated traffic conditions, 
channel error propagation and capture effects over a Rayleigh 
fading channel under the hypothesis of employing a 2-way 
handshaking access mechanism. 

It is useful to briefly recall the 2-way handshaking mech- 
anism. A station that wants to transmit a packet, waits until 
the channel is sensed idle for a Distributed InterFrame Space 
(DIES), follows the backoff rules and then transmits a packet. 
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Fig. 1. Markov chain for tlie contention model in unsaturated traiSc 
conditions, based on the 2-way handshaking technique, considering the effects 
of capture and channel induced errors. 

After the successful reception of a data frame, the receiver 
sends an ACK frame to the transmitter Only upon a correct 
ACK frame reception, the transmitter assumes successful 
delivery of the corresponding data frame. If an ACK frame is 
received in error or if no ACK frame is received, due possibly 
to an erroneous reception of the preceding data frame, the 
transmitter will contend again for the medium. 

On the basis of this assumption, collisions can occur with 
probability Pcoi on the transmitted packets, while transmission 
errors due to the channel, can occur with probability P^. 
We assume that collisions and transmission error events are 
statistically independent. In this scenario, a packet is suc- 
cessfully transmitted if there is no collision (this event has 
probability 1 — Pcoi) and the packet encounters no channel 
errors during transmission (this event has probability 1 — Pg). 
The probability of successful transmission is therefore equal 
to (1 — Pe)(l — Pcoi), from which we can set an equivalent 
probability of failed transmission as P^q = Pe + Pcoi — PePcoi- 

Furthermore, in mobile radio environment, it may happen 
that the channel is captured by a station whose power level is 
stronger than other stations transmitting at the same time. This 
may be due to relative distances and/or channel conditions for 
each user and may happen whether or not the terminals exer- 
cise power control. Capture effect often reduces the collision 
probability on the channel since the stations whose power level 
at the receiver are very low due to path attenuation, shadowing 
and fading, are considered as interferers at the access point 
(AP) raising the noise floor 

To simplify the analysis, we make the assumption that the 
impact of the channel induced errors on the packet headers 
are negligible because of their short length with respect to 
the data payload size. This is justified on the basis of the 
assumption that the bit errors affecting the transmitted data 



are independent of each other. Hence, the packet or frame 
error rate, identified respectively with the acronyms PER or 
FER, is a function of the packet length, with shorter packets 
having exponentially smaller probability of error compared to 
longer packets. 

Practical networks operate in unsaturated traffic conditions. 
In this case, Bianchi's model ifTSl assuming the presence of a 
packet to be transmitted in each and every station's buffer, is 
not valid anymore. However, the simplicity of such a model 
can be retained also in unsaturated conditions by introducing 
a new state, labelled /, accounting for the following two 
situations: 



« immediately after a successful transmission, the buffer of 

the transmitting station is empty; 
• the station is in an idle state with an empty buffer until 

a new packet arrives at the buffer for transmission. 



With these considerations in mind, let us discuss the Markov 
model shown in Fig.[Tl modelling unsaturated traffic condition. 
Notice that in our Markov chain we do not model the 802.11 
post-backoff feature. However, post-backoff has negligible 
effects on the theoretical results developed in the paper as 
exemplified by the simulation results obtained in the following. 
Similar to the model in [13], (m+ 1) different backoff stages 
are considered (this includes the zero-th stage). The maximum 
contention windows (CW) size is CWmax = 2"^W, and 
the notation Wi — is used to define the z*'* contention 
window size. A packet transmission is attempted only in the 
(i, 0) states, Vi = 0, . . . , m. If collision occurs, or transmission 
is unsuccessful due to channel errors, the backoff stage is 
incremented, so that the new state can be (i + 1,/c) with 
probability Peq/VFi+i, since a uniform distribution between 
the states in the same backoff stage is assumed. We consider 
capture as a subset of the event of a collision (i.e., capture 
implicitly implies a collision). In other word, the capture 
event can happen in the presence of collision by allowing the 
transmitting station with the highest received power level at the 
access point to capture the channel. In the event of capture, 
no collision is detected and the Markov model transits into 
one of the transmitting states (i,0) depending on the current 
contention stage. If no collision occurs (or is detected), a data 
frame can be transmitted and the transmitting station enters 
state (i, 0) based on its backoff stage. From state (i, 0) the 
transmitting station re-enters the initial backoff stage i = if 
the transmission is successful and at least one packet is present 
in the buffer, otherwise the station transit in the state labelled 
/ waiting for a new packet arrival. Otherwise, if errors occur 
during transmission, the ACK packet is not sent, an ACK- 
timeout occurs, and the backoff stage is changed to (i + 1, fc) 
with probability Pgg/W'i+i- 

The Markov Process of Fig. [T] is governed by the following 



transition probabilities 



= 1, 



fee [o,w^-2], ie [o,m] 



Po,k\^,o = 9(1 - ^'e5)/W^o, k e [0, Wo-1], ie [0,m] 
PiMi-T.fi =Peq/Wi, fc G [0,W^, - 1], i G [l,m] 

Pn^.k\m.O = P.JW^, k G [0, - 1] 

f'n^.o = (l-g)(l-Pe,), iG[0,m] 

^'o.fcl/ =g/W^o, fcG[0,M/o-l] 
= 1 - ? 

(1) 

The first equation in ^ states that, at the beginning of each 
slot time, the backoff time is decremented. The second equa- 
tion accounts for the fact that after a successful transmission, 
a new packet transmission starts with backoff stage with 
probability q, in case there is a new packet in the buffer to be 
transmitted. Third and fourth equations deal with unsuccessful 
transmissions and the need to reschedule a new contention 
stage. The fifth equation deals with the practical situation 
in which after a successful transmission, the buffer of the 
station is empty, and as a consequence, the station transits 
in the idle state / waiting for a new packet arrival. The sixth 
equation models the situation in which a new packet arrives in 
the station buffer, and a new backoff procedure is scheduled. 
Finally, the seventh equation models the situation in which 
there are no packets to be transmitted and the station is in the 
idle state. 

III. Markovian Process Analysis and Throughput 
Computation 

Next line of pursuit consists in finding a solution of the 
stationary distribution: 

h^k = lim P[s{t) = i,h{t) = fc], Vfc G [0, VK,-1], Vi G [0, m] 

that is, the probability of a station occupying a given state at 
any discrete time slot. 

First, we note the following relations: 

bifi = Peq ■ bi-ifl • 6o,o, Vi G [1, m - 1] 

Om,0 = i_p^^ • oo,o, I = m 

whereby, P^q El is the equivalent probability of failed transmis- 
sion, that takes into account the need for a new contention due 
to either packet collision (Pcoi) or channel errors (Pe), i.e., 

Peq = Pool +Pe-Pe- Pool 0) 

State bi in Fig. [T] considers both the situation in which after a 
successful transmission there are no packets to be transmitted, 
and the situation in which the packet queue is empty and 
the station is waiting for new packet arrival. The stationary 
probability to be in state bj can be evaluated as follows: 

hi - {I - q)ii- Peq) ■Etoho + a -q)-bi 

_ (l-g)(l-Pe,) , (4) 

- q • Z^i=0 "i.O 

The expression above reflects the fact that state bj can be 
reached after a successful packet transmission from any state 

'■Pi,fe|j.n is short for P{s(t + 1) = i, b{t + 1) = k\s{t) = j, b{t) = n}. 

^For simplicity, we assume that at each transmission attempt any station 
wiU encounter a constant and independent probabiHty of failed transmission, 
Peq, independently from the number of retransmissions already suffered from 
each station. 



bifl, Vi G [0, m] with probability (1 — q){l — Peq), or because 
the station is waiting in idle state with probability (1 — q), 
whereby q is the probability of having at least one packet to 
be transmitted in the buffer. The statistical model of q will be 
discussed in the next section. 

The other stationary probabilities for any k G — 1] 

follow by resorting to the state transition diagram shown in 

Fig.m 



q{l - Peq) ■ X^Ho ^'^0 +<l-bl, 
Peq • &i-l,0 5 
Peq{bm-lfi + &m,o): 



i = 

i G [1, m — 1] 
i — m 

(5) 



Upon substituting dU in (|5]l, 6o,fc can be rewritten as follows: 

q{l - Peq) ■ E"o 

= - Peq) ■ TZo ho + q ■ ii^^Hi^ . 6,^ „ 

= {l-Peq)[q+{l-q)]-j:toho 

= i^-Peq)■ET=ob^M 

(6) 

Employing the normalization condition, after some mathemat- 
ical manipulations, and remembering the relation X]"=o ^*.o = 



bo.o 
1-P„ 



it is possible to obtain: 



1 = T,T=oT,Y^o^Kk + bi 



+ 



(1-9)(1-Peg) 



(l-q)(l-Peg) 6o,o 
^ q 



2Pe,)" 



2(1-9) 
9 

a 

Normalization condition yields the following equation for 
computation of bo o: 
^ 2 

^ 2(1-Pe,) 

_ 2(l-Pe,)(l-2P,,)g 

- q[(W+l)(l-2P,,)+WPe,(l-(2P,,)")]+2(l-g)(l-Pe,)(l-2Pe,) 

(8) 

As a side note, if q ^ 1, i.e., the station is approaching 
saturated traffic conditions and assuming packet transmission 
errors are only due to collisions, i.e., Pe = and Peq = Pcoi, 
dHJ we get: 

2(l-Pcoi)(l-2P,o,) 



U . 2(l-P,oi)(l-2P,oi) ^QS 

imig^l uo,0 ^ (Vi/+i)(i_2P^„,)+WPeoi(l-(2fcoi)'") 

which corresponds to the stationary state probability 5o,o found 
by Bianchi ifTSl under saturated conditions. 

Equ. ^ is then used to compute t, the probability that a 
station starts a transmission in a randomly chosen time slot. 
In fact, taking into account that a packet transmission occurs 
when the backoff counter reaches zero, we have: 

^ — v-^m , _ bo.o 

„ 2(l-2Pe,)g 

~ <;l(W+l)(l-2Pe,)+WP,,(l-(2P,,)'")l+2(l-g)(l-Pe,)(l-2Pe,) 

(10) 

Once again as a consistency check, note that if to = 0, 
i.e., no exponential backoff is employed, g — > 1, i.e., the 
station is approaching saturated traffic conditions, and under 



the condition that packet transmission errors are only due to 
colUsions, i.e., Pg = and P^q = Pcoi, ( flOl t we get: 



lim T 



W+1 



(11) 



which is the result found in 11261 for the constant backoff 
setting, showing that the transmission probability r is inde- 
pendent of the collision probability Pcoi- 

The collision probability needed to compute r can be found 
considering that using a 2-way hand-shaking mechanism, a 
packet from a transmitting station encounters a collision if 
in a given time slot, at least one of the remaining {N ~ 1) 
stations transmits simultaneously another packet, and there is 
no capture. In our model, we assume that capture is a subset 
of the collision events. This is indeed justified by the fact 
that there is no capture without collision, and that capture 
occurs only because of collisions between a certain number 
of transmitting stations attempting to transmit simultaneously 
on the channel. 



Pcoi = 1 - (1 - r) 



-P 



cap 



(12) 



As far as the capture effects are concerned, we resort to the 
mathematical formulation proposed in ll23l . Il24l . 

We consider a scenario in which N stations are uniformly 
distributed in a circular area of radius R, and transmit to- 
ward a common access point placed in the center of such 
an area. When signal transmission is affected by Rayleigh 
fading, the instantaneous power of the signal received by the 
receiver placed at a mutual distance from the transmitter is 
exponentially distributed as: 

1 

f{x) — — e po , x > 

Po 

whereby po is the local-mean power determined by the equa- 
tion: 

Po^A- r'T'^'-Pt^ (13) 

where, Up is the path-loss exponent (which is typically greater 
or equal to 3.5 in indoor propagation conditions in the absence 
of the direct signal path), Ptx is the transmitted power, and 
A ■ r^"'' is the deterministic path-loss |29] . Both A and Pt^ 
are identical for all transmitted frames. 

Under the hypothesis of power-controlled stations in infras- 
tructure mode and Rayleigh fading, the capture probability 
conditioned on i interfering frames can be defined as follows: 



Pep (7 > Zog{Sf )\i) 



1 



[1 + ^o9{Sf)T 



whereby, 7, defined as 



(14) 



(15) 



is the ratio of the power of the useful signal and the sum of 
the powers of the i interfering channel contenders transmitting 
simultaneously i frames, g{Sf) is the inverse of the processing 
gain of the correlation receiver, and zq is the capture ratio, i.e., 
the value of the signal-to-interference power ratio identifying 
the capture threshold at the receiver Notice that ( fT4] i signifies 



the fact that capture probability corresponds to the probability 
that the power ratio 7 is above the capture threshold zog{Sf) 
which considers the inverse of the processing gain g{Sf). For 
Direct Sequence Spread Spectrum (DSSS) using a 11-chip 
spreading factor (5/ = 11), we have g{Sf) = EM . 

Upon defining the probability of generating exactly i + 1 
interfering frames over N contending stations in a generic 
slot time: 

^ V'^+i(l_r)^-'-i 



vi + 1, 

the frame capture probability Pcap can be obtained as follows: 



P. 



N-l 

cap — 

i— 



^ G + 1)^^^'^^ rf-'Pcp (7 > ZogiSfm 



(16) 

Putting together Equ.s (O, (fTOl l. (1% . and (fT6] l the following 
nonlinear system can be defined and solved numerically, 
obtaining the values of r, Pcoi, Pcap, and Pg^: 

^ _ 2(l-2P,,)g 

g[(VK+l)(l-2Pe„) + VKPe<,(l-(2P,,)-)]+2(l-g)(l-P,,)(l-2Pe,) 

Pcoi = 1 - (1 - r)^-l ~ Pcap 
Peq Pcoi Pc Pe ' Pcoi 



Pcap-EtiL^,y+Hi-ry 



(l+zogiSf))' 



(17) 

The final step in the analysis is the computation of the 
normalized system throughput, defined as the fraction of time 
the channel is used to successfully transmit payload bits: 

Pt ■ P. ■ (1 - Pe)E[PL] 

(1 - Pt)a + Pt(l - Ps)T, + PPs(l - Pe)T, + PtPsPeT, 

(18) 

where, 

• Pt is the probability that there is at least one transmission 
in the considered time slot, with N stations contending 
for the channel, each transmitting with probability r: 



Pt = 1 - (1 - r) 



N 



(19) 



Ps is the conditional probability that a packet transmis- 
sion occurring on the channel is successful. This event 
corresponds to the case in which exactly one station 
transmits in a given time slot, or two or more stations 
transmit simultaneously and capture by the desired station 
occurs. These conditions yields the following probability: 



P. 



Pc, 



Pt 



(20) 



Tc, Te and P, are the average times a channel is 
sensed busy due to a collision, error affected data frame 
transmission time and successful data frame transmission 
times, respectively. Knowing the time durations for ACK 
frames, ACK timeout, DIFS, SIFS, a, data packet length 
(PL) and PHY and MAC headers duration (H), and 
propagation delay Tp, Tc, Tg, and Tc can be computed 
as follows |17|: 

Tc = H + PL + AC Ktimeout 
Tc = H + PL + AC Ktimeout 

T, =H + PL + SIFS + Tp + ACK + DIFS + Tp 

(21) 



> E[PL] is the average packet payload length. 
• fj is the duration of an empty time slot. 
The setup described above is used in Section |IV] for DCF 
simulation at the MAC layer. 

A. Modelling offered load and estimation of probability q 

In our analysis, the offered load related to each station is 
characterized by parameter A representing the rate at which 
packets arrive at the station buffer from the upper layers, and 
measured in pkt/s. The time between two packet arrivals is 
defined as interarrival time, and its mean value is evaluated 
as J. One of the most commonly used traffic models assumes 
packet arrival process is Poisson. The resulting interarrival 
times are exponentially distributed. 

In the proposed model, we need a probabiUty q that indicates 
if there is at least one packet to be transmitted in the queue. 
Probability q can be well approximated in a situation with 
small buffer size 1271 through the following relation: 

q = l- e^-^-^I^'^l (22) 

where, E[Sts] is the expected time per slot, which is useful to 
relate the state of the Markov chain with the actual time spent 
in each state. 

A more accurate model can be derived upon considering 
different values of q for each backoff state. However, a rea- 
sonable solution consists in using a mean probability valid for 
the whole Markov model |27|, derived from E[Sts]- The value 
of E[Sts] can be obtained by resorting to the durations and 
the respective probabilities of the idle slot (a), the successful 
transmission slot (Tg), the error slot due to collision (Tc), and 
the error slot due to channel (Tg), as follows: 

E[St,] - il-Pt)-a + Ptil-Ps)-T,+ 

+PtPsPe ■ T, + PtPsil -Pe)-Ts 

Upon recalling that packet inter-arrival times are exponentially 
distributed, we can use the average slot time to calculate 
the probability q that in such a time interval a given station 
receives a packet from upper layers in its transmission queue. 
The probability that in a generic time T, k events occur, is: 

P{aiT) = k} = e-'^^ (24) 
from which we obtain the relation (l22ll referred to earlier: 

q = l- P{a{E[Sts]) = 0} = 1 - e-^^[^-l (25) 

IV. Simulation Results and Model Validations 

This section focuses on simulation results for validating the 
theoretical models and derivations presented in the previous 
sections. We have developed a C++ simulator modelling 
both the DCF protocol details in 802.11b and the backoff 
procedures of a specific number of independent transmitting 
stations. The simulator is designed to implement the main 
tasks accomplished at both MAC and PHY layers of a wireless 
network in a more versatile and customizable manner than ns- 
2, where the lack of a complete physical layer makes difficult 
a precise configuration at this level. 



TABLE I 
Typical network parameters 



MAC header 


24 bytes 


PHY header 


16 bytes 


Payload size 


1024 bytes 


ACK 


14 bytes 


RTS 


20 bytes 


CTS 


14 bytes 


Tp 


1 flS 


Slot time 


20 /IS 


SIFS 


10 lis 


DIPS 


50 /IS 


EIFS 


300 fis 


ACK timeout 


300 fis 


CTS timeout 


300 fis 



The simulator considers an Infrastructure BSS (Basic Ser- 
vice Set) with an AP and a certain number of mobile stations 
which communicates only with the AP. For the sake of 
simplicity, inside each station there are only three fundamental 
working levels: traffic model generator, MAC and PHY layers. 
Traffic is generated following the exponential distribution for 
the packet interarrival times. Moreover, the MAC layer is 
managed by a state machine which follows the main directives 
specified in the standard [11, namely waiting times (DIFS, 
SIFS, EIFS), post-backoff, backoff, basic and RTS/CTS access 
mode. 

Typical MAC layer parameters for IEEE802.11b are given 
in Table U [ 1] . In so far as the computation of the FER is 
concerned, it should be noted that data transmission rate of 
various packet types differ. For simplicity, we assume that data 
packets transmitted by different stations are affected by the 
same FER. This way, channel errors on the transmitted packets 
can be accounted for as it is done within ns-2 |30|. In other 
words, a uniformly distributed binary random variable is 
generated in order to decide if a transmitted packet is received 
erroneously. The statistic of such a random variable is P{Xe = 
1) = Pe{SNR) (as specified in dH), and P{Xe = 0) = 
1-PeiSNR). 

As far as capture effects are concerned, the investigated 
scenario is as follows. In our simulator, N stations are ran- 
domly placed in a circular area of radius R (in the simulation 
results presented below we assume R = 10m), while the 
AP is placed at the center of the circle. When two or more 
station transmissions collide, the value of 7 as defined in (flST l 
is evaluated for any transmitting station given their relative 
distance r,; from the AP. Let 7^ be the value of 7 for the 
j-th transmitting station among the i + 1 colliding stations. 
The power between a transmitter and a receiver is evaluated in 
accordance to (fTsT l with a path-loss exponent equal to 3.5. The 
values of jj for each colliding station are compared with the 
threshold Zq ■ g{Sf). Then, the transmitting station for which 
jj is above the threshold captures the channel. 

The physical layer (PHY) of the basic 802.11b standard is 
based on the spread spectrum technology. Two options are 
specified, the Frequency Hopped Spread Spectrum (FHSS) 
and the Direct Sequence Spread Spectrum (DSSS). The FHSS 
uses Frequency Shift Keying (FSK) while the DSSS uses 



Differential Phase Shift Keying (DPSK) or Complementary 
Code Keying (CCK). The 802.11b employs DSSS at vaiious 
rates including one employing CCK encoding 4 and 8 bits on 
one CCK symbol. The four supported data rates in 802.11b 
are 1, 2, 5.5 and 11 Mbps. 

In particular, the following data are transmitted at lowest 
rate (1 Mbps) in IEEE802.11b: PLCP=16 bytes (PLCP plus 
header), ACK Headers=16 bytes, ACK=14 bytes, RTS=20 
bytes, CTS=14 bytes (RTS and CTS are only for four way 
handshake). 

The PER as a function of the SNR can be computed as 
follows: 

P^SNR) = 1- [1 - P4PLCP, SNR)] ■ [1 - P,{DATA, SNR)] 

(26) 

where, 

Pe{PLCP, SNR) = 1 - [1 - Pb{BPSK, SNR)]^""^^^^ , (27) 
and 

P,{DATA,SNR) = - Pb{TYPE,SNR)]^''^'^'^'^^+''''^'^^ . 

(28) 

Pb{BPSK,SNR) is the BER as a function of SNR for 
the lowest data transmit rate employing DBPSK modulation, 
DATA denotes the data block size in bytes, and any other 
constant byte size in above expression represents overhead. 
Note that the FER, Pe{SNR), implicidy depends on the 
modulation format used. Hence, for each supported rate, one 
curve for Pe{SNR) as a function of SNR can be gener- 
ated. Pb{TYPE, SNR) is modulation dependent whereby the 
parameter TYPE can be any of the following TYPE G 
{DBPSK, DQPSK, CCK5.5, CCKllj^ 

For DBPSK and DQPSK modulation formats, 
Pb{TYPE,SNR) can be well approximated by ll28l : 

max(^,l) 

2 V i r - de 

max(log2M,2) ^ n i + sin^ ( ^2l^) 

(29) 

whereby, M is the number of bits per modulated symbols, 7 
is the signal-to-noise ratio, and 6 is the signal direction over 
the Rayleigh fading channel. 

We note that the proposed DCF model is valid with any 
other PHY setup. Actually, all we need is simply the packet 
error rate probabiHty Pe{DATA,SNR) in ^ for the spe- 
cific PHY and channel transmission conditions model. All the 
mathematical derivations proposed in this paper are specified 
with respect to Pe{DATA, SNR), and so they are vaUd for 
any kind of transmission model at the PHY layer. 

In what follows, we shall present theoretical and simulation 
results for the lowest supported data rate. We note that by 
repeating the process, similar curves can be generated for 
all other types of modulation formats. All we need is really 
the BER as a function of SNR for each modulation format 
and the corresponding raw data rate over the channel. If the 
terminals use rate adaptation, then under optimal operating 
condition, the achievable throughput for a given SNR is the 

'The acronyms are short for Differential Binary Phase Shift Keying, 
Differential Quadrature Phase Shift Keying and Complementary Code Keying, 
respectively. 
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Fig. 2. Theoretical and simulated throughput for the 2-way mechanism as 
a function of the SNR, for two different packet rates and capture thresholds, 
shown in the legend. Simulated points are identified by star-markers over the 
respective theoretical curves. Payload size is 1024 bytes, while the number of 
contending stations is Af = 10. 
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Fig. 3. Theoretical and simulated throughput for the 2-way mechanism as a 
function of the packet rate A, for three different number of contending stations 
and two different values of SNR as noted in the legend. Capture thi'esholds is 
zq = 6dB. Simulated points are identified by star-markers over the respective 
theoretical curves. Payload size is 1024 bytes. 



maximum over the set of modulation formats supported. We 
have verified a close match between theoretical and simulated 
performances for other transmitting data rates as well. In 
the results presented below we assume the following values 
for the contention window: CWmin = 32, m = 5, and 

CWrna. = 2™ • CW,mn = 1024. 

Fig. 12] shows the behavior of the throughput for the 2- 
way mechanism as a function of the SNR, for two different 
capture thresholds, zq, and packet rates. A, as noted in the 
legend, and for = 10 transmitting stations. Simulated 
points are marked by star on the respective theoretical curves. 
Throughput increases as a function of the three parameters 
SNR, A, and zq until saturation point. 

Fig. [3] shows the behavior of the throughput as a function 
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Fig. 4. Theoretical and simulated throughput for the 2-way mechanism as 
a function of the packet rate A, for three different number of contending 
stations and two different values of SNR as noted in the legend. Capture 
threshold is zq = 24dB. Simulated points are identified by star-markers over 
the respective theoretical curves. Payload size is 1024 bytes. 
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Fig. 5. Theoretical and simulated throughput for the 2-way mechanism as a 
function of the packet rate A, for three different number of contending stations 
and two different capture thresholds noted in the legend. Simulated points are 
identified by star-markers over the respective theoretical curves. Payload size 
is 128 bytes, while SNR=45dB. 



of A, i.e., the packet rate, for three different values of the 
number of contending stations and for two values of SNR. The 
capture threshold is zq = 6dB. Beside noting the throughput 
improvement achievable for high SNR, notice that for a spec- 
ified number of contending stations, the throughput manifests 
a linear behavior for low values of packet rates with a slope 
depending mainly on the number of stations N. However, for 
increasing values of A, the saturation behavior occurs quite 
fast. Notice that, as exemplified in (IZST l. g ^ 1 as A oo. 
Actually, saturated traffic conditions are achieved quite fast 
for values of A on the order of ten packets per second with a 
number of contending stations greater than or equal to 10. 

Fig. |4] shows the behavior of the throughput as a function 
of A, for three different values of the number of contending 
stations and for two different SNRs. The capture threshold 
is Zq = 24dB. We can draw conclusions similar to the ones 
derived for Fig.|3] Upon comparing the curves shown in Figs [3] 
and |4] it is easily seen that capture effects allow the system 
throughput to be almost the same independently from the 
number of stations in saturated conditions, i.e., for high values 
of A. 

Fig Elalso shows the presence of a peak in the throughput as 
a function of A, which characterizes the transition between the 
linear and saturated throughput. Such a peak tends to manifest 
itself for increasing values of A as the number of stations 
N increases. A comparative analysis of the curves shown in 
Figs [3] and |4] reveals that the peak of the throughput tends to 
disappear because of the presence of capture effects during 
transmission. 

Notice also that the throughput behavior exhibits saturation 
quite fast for low values of SNR independently from capture 
effects; this is essentially due to the fact that channel propaga- 
tion errors tend to dominate over both collisions and capture 
when the SNR is low. 

In order to assess throughput performances as a function of 
the payload size. Fig. |5] shows the behavior of the throughput 



as a function of A, for a payload size equal to 128 bytes. The 
other simulation parameters are noted in the legend. Beside 
noting similar conclusions on the effects of both capture and 
number of contending stations over the throughput behavior as 
for the previous mentioned figures, notice also that a smaller 
payload size reduces proportionally the slope of the throughput 
before saturated conditions, and the maximum achievable 
throughput in saturated conditions. This issue will be dealt 
with in the next section. 

A. Simulation results within ns-2 

For the sake of validating the proposed model and our 
developed simulator against a widely adopted network sim- 
ulator, while providing simulation results for a new set of 
parameters, we have performed new simulations using ns-2 
(version 2.29) for throughput evaluations. Due to the lack 
of a complete PHY layer implementation, we made some 
modifications in the original ns-2 source code in order to 
account for a specific FER on the transmitted packets. In 
this respect, we adopted the suggestions proposed in [30| for 
simulating a 802.11b channel within ns-2. Briefly, channel 
errors on the transmitted packets are accounted for by using 
a uniformly distributed binary random variable with the 
following probabilities: P{Xe = 1) = Pe{SNR) (as specified 
in and P{Xe = 0) = 1 - P^{SNR). For simpHcity, we 

assume that data packets transmitted by different stations are 
affected by the same FER. 

The employed propagation model is identified by the label 
Propagation/Shadowing within ns-2. We considered a 
path-loss exponent equal to 3.5 with a zero standard deviation 
in compliance with the theoretical capture scenario described 
above. Notice that the standard ns-2 simulator provides only 
Ad-Hoc support and related routing protocols; for this reason 
we applied a patch for implementing an Infrastructure BSS 
without any routing protocol between the various involved 
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Fig. 6. Theoretical and simulated (ns-2) throughput for the 2-way mechanism 
as a function of the packet rate A, for two different payload sizes and packet 
error probabilities, shown in the legend, and without capture. Simulated points 
are identified by cross-markers over the respective theoretical curves. The 
number of contending stations is = 10. 

Stations. This patch is identified by the acronym NOAH, 
standing for NO Ad-Hoc, and can be downloaded from 
http://icapeople.epfl.ch/widmer/uwb/ns-2/noah/. 

In connection with the employed traffic model, we used 
an exponential distribution for simulating packet interarrival 
times. 

Throughput evaluation is accomplished by averaging over 
100 sample scenarios, whereby any transmitting scenario 
considers a set of N randomly distributed (with a uniform 
pdf) stations over a circular area of radius 10 m. Simulation 
results along with theoretical curves are shown in Figs |6][8] 
The employed parameters are noted in the respective legends 
of each figure. While the main conclusions related to the 
throughput behavior as a function of the number of contending 
stations, capture threshold, traffic load, and FER, are similar 
to the ones already derived above with our C++ simulator, 
here we notice a good agreement between theoretical and ns- 
2 simulation results for a wide range of packet arrival rates, 
A. 

Figs |6]|7] also depict the theoretical saturated throughput 
derived by Bianchi in llT3l in order to underline a good 
matching between our theoretical model, valid for any traffic 
load in the investigated range of A up to the saturation 
condition for A — > oo. 

V. Modelling the linear behavior of the 

THROUGHPUT IN UNSATURATED CONDITIONS 

The linear behavior of the throughput in unsaturated traffic 
conditions, along with its dependence on some key network 
parameters, can be better understood by analyzing the through- 
put in ( fTSl l as A ^ 0, i.e., in unloaded traffic conditions. 
From ( fTOl ). it is easily seen that the probability t ^ when 
A ^ 0, since g ^ 0. In particular, for very small values of A, 
the following approximation from (fTOl l holds: 

q 

T ~ 

l-Pe 
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Fig. 7. Theoretical and simulated (ns-2) throughput for the 2-way mechanism 
as a function of the packet rate A, for two different payload sizes and packet 
error probabilities, shown in the legend, and without capture. Simulated points 
are identified by cross-markers over the respective theoretical curves. The 
number of contending stations is A'^ = 5. 
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Fig. 8. Both subplots show theoretical and simulated (ns-2) throughput for 
the 2-way mechanism as a function of the packet rate A, for two different 
payload sizes and fixed packet error probabihty, as noted in the legend. 
Capture threshold zq = IdB in both subplots. Simulated points are identified 
by cross-markers over the respective theoretical curves. Upper subplot refers 
to a number of contending stations equal to N = 10, whereas the lower 
subplot corresponds to A'^ = 5. 

In addition, it is straightforward to verify that Pcap 
(from (fT6t). As A approaches zero, the throughput can be 
approximated by: 

a a 

(30) 

since q = \ — ^-^E\Sts\ ^.^jj approximated by g ~ E\Sts\X 
upon employing MacLaurin approximation of the exponential, 
and i?[S'ts] cr for A ^ as it can be deduced from (|23]) . 
The main conclusion that can be drawn from this equation 
is that as A — > 0, the throughput tends to manifest a linear 
behavior relative to the packet rate A with a slope depending 
on both the number of contending stations and the average 



Values of Ac from qI 



TABLE II 

I FOR THE SETUP CONSIDERED IN FIGS|4I5I 



E[PL]=I024 byte 


N=4 


N=10 


N=20 


SNR=20dB 


13.9187 


5.5372 


2.7638 


SNR=45dB 


26.7546 


10.6444 


5.3132 


E[PL]=128 byte 


N=4 


N=10 


N=20 


SNR=45dB 


135.0307 


53.3990 


26.6039 



pay load size E[PL]. 

It is interesting to estimate the interval of validity [0, Ac] of 
the linear throughput model proposed in ( |30] |. To this end, let 
us rewrite ( fTSl ) as follows: 



S = 



E\PL] 



(31) 



1 - Pe ^ 1 - Pe 



Ps(l-Pc) 

Since Ts,Tc,Te and Pg are independent of r, S can be 
maximized by minimizing the function: 

Fir) = — (32) 



Psil-Pe) 

or equivalently, by maximizing the function: 

1 _ Ps{l-Pe) 



Flir) = 



F{t) 



- i-Pt 
Pt 



Tr 



(33) 



with respect to r. Upon substituting (fT9] l and ( |20| | in Fi{t) 
it is possible to write: 

Nt{1 - Pe) 



(34) 



fT(l - r) + rc(l - - Tc(l - r) 

Differentiating Fi (r) with respect to r and equating it to zero 
it is possible to write: 



{l-Tr{a~T,)~T,{NT-l) = 



(35) 



Under the hypothesis r <C 1, the following approximation 
holds: 



(1 



1-Nt + 



N{N ~l) 



air) 



By substituting the previous approximation in ( |35] |. it is 
possible to obtain the value of t maximizing Fi{t): 



_ a-^a [Na-2{N-l)ia-T,)]/N 

Substituting t„i in dSTT i. it is possible to obtain the maximum 
throughput Sm- 

E[PL] 



Ts - 



1-Pe 



((T-Tj(l-T„)"+Te 

Afr,„(l-r„)""i(l-Pe) 



(37) 



Equating S*™ to the linear model in (l30l l. it is possible to obtain 
the value Ac for which the linear model reaches the beginning 
of the saturation zone, i.e., the maximum A above which the 
linear model loses validity: 



Ac 



(38) 



The linear model dSOl l shows a close agreement with both 
theoretical and simulated throughput curves. Fig. |9] shows the 
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Fig. 9. Throughput for the 2-way mechanism as a function of the packet 
rate A, for thi'ee different number of contending stations. Capture threshold 
and SNR are set to very high values since in pre-saturation regime. Payload 
size is 102 4 by tes. Straight lines refer to the linear model of the throughput 
derived in 130). 



Straight lines (l30l l for three different values of the number of 
contending stations N, compared with the theoretical curves 
akeady depicted in Fig. |4] The figure also shows the values 
of S„i as derived from (IJTT i along with the three values of Ac 
deduced from (l38T l. Notice that Ac, besides being the limit of 
validity of the linear model, can also be interpreted as the value 
of packet rate above which the throughput enters saturation 
conditions. 

In order to further verify the values of Ac deduced from ( |38] ). 
Table HIl shows the values of Ac related to both simulation and 
theoretical results depicted in Figs |4]|5] In particular, values 
shown in the upper part of the table, i.e., the ones labelled 
E[PL]=1024, are related to the curves in Fig. |4] whereas the 
other values are related to the results in Fig. |5] Notice the 
good agreement between theoretical and simulated results even 
in the presence of capture. The key observation here can be 
deduced by a comparative analysis of the results shown in 
Figs. [3]|4| capture tends to simply shape the peak on which 
the throughput reaches its maximum, leaving the abscissa of 
the maximum, that is Ac, in the same position. 

VI. Conclusions 

In this paper, we have provided an extension of the Markov 
model characterizing the DCF behavior at the MAC layer of 
the IEEE802.il series of standards by accounting for channel 
induced errors and capture effects typical of fading envi- 
ronments under unsaturated traffic conditions. The modelling 
allows taking into consideration the impact of channel con- 
tention in throughput analysis which is often not considered 
or it is considered in a static mode by using a mean contention 
period. Subsequently, based on justifiable assumptions, the 
stationary probability of the Markov chain is calculated to 
obtain the behavior of the throughput in both unsaturated 
and saturated conditions. The closed form expressions allow 
derivation of the throughput as a function of a multitude of 
system level parameters including packet and header sizes for a 



variety of applications. Simulation results confirm the validity 
of the proposed theoretical models. 
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